DESY 95-241 
December 1995 



On the Theoretical Status of 
Deep Inelastic Scattering 



J. Bliimlein 



DESY-Zeuthen, 
Platanenallee 6, D-15735 Zeuthen, Germany 



Abstract 

The theoretical status of perturbative QED and QCD corrections to deep inelastic scatter- 
ing is reviewed. 



Invited talk at the XXV International Symposium on Multiparticle Dynamics, 
Stara Lesna, Slovakia, 11-17 September, 1995 



ON THE THEORETICAL STATUS OF 
DEEP INELASTIC SCATTERING 



JOHANNES BLUMLEIN 
DESY-Zeuthen, Platanenallee 6, 

D-15735 Zeuthen, Germany 
e-mail: HlKBLU@dsyibm.desy.de 



ABSTRACT 

The theoretical status of perturbative QED and QCD corrections to deep inelas- 
tic scattering is reviewed. 



1. Introduction 

2. QED radiative corrections to DIS 

3. The running coupling constant 

4. The evolution equation 

4.1. Splitting functions 

4.2. Coefficient functions 

4.3. 0{a^g) corrections 

5. Resummation of small x contributions 

5.1. Singlet terms 

5.2. Non-singlet terms 

6. Heavy flavour contributions to structure functions 

7. J/^ production 

8. QCD corrections to polarized structure functions 

9. Open problems 



1. Introduction 



Deep inelastic lepton-hadron scattering provides one of the cleanest ways to investi- 
gate the nucleon structure at short distances. Since both charged leptons (e^,/i^) and 
neutrinos {y-fV) may be used as probes in neutral and charged current deep inelastic 
scattering off protons or isoscalar targets a variety of scattering cross sections can be 
measured containing different flavour combinations of quarks. Furthermore, the spin 
structure of nucleons can be investigated using both polarized leptons and targets. 

The basic diagram describing the process is shown in figure 1. The neutral and 
charged current scattering cross sections are given by 

f^_p . 2^ MNSy B 

in the Born approximation. The leptonic and hadronic tensors L'^,W^^ depend 
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on the quantum numbers of the lepton and hadron, respectively, and those of the 
exchanged boson ^. Here, x and y denote the Bjorken variables, = {I — Z')^, 
and s = {I -\- Py. The factor Pi^b{Q^) collects couplings and propagator terms, e.g. 
Pi±,-,{Q^) = 2TvayQ\ P,,w+{Q^) = GIM^/{Atv{Q^ + M^f), with a being the fine 
structure constant, Gp the Fermi constant, and Mw the W boson mass. 




Figure 1: Diagram describing deep inelastic IN scattering, with B = f,Z,W^, N = p,n, and 
/ = l^,v,V. 

The hadronic tensor may be represented by the structure functions describing the 
scattering process. In the parton model the structure functions are expressed in 
terms of parton densities, qi[x,Q^). For l^N scattering the structure functions are 
thus given by: 



xGs{x,Q') 
xHs{x,Q') 
W+{x,Q') 
W,-{x,Q') 
xW+{x,Q') 
xWi{x,Q') 



xY^el[q{x,Q') + q{x,Q')] \^\' 

2xY,e,vM^,Q')+t^,Q')] 

xY.{vl + alMx,Q') + q{x,Q')] \Z\^ 

xYvqa^\q{x,Q'^)-q(x,Q'^)\ \Z\^ 

2xY\Ai{x,Q-')\ui{x,Q-')\ \W^\^ 

2xY\ni{x,Q-')\2i{x,Q-')\ \W-\^ 

1xY\i,{x,Q-')-u,{x,Q-')\ \W^\^ 

2xYXni{x,Q-')-2i{x,Q-')\ \W-\^ 



(2) 



where eq,Vq, and denote the charge, vector-, and axialvector coupling constants 
of the quarks. The u[u)N scattering processes are described by six further structure 
functions W^'^ , xW^'^ , F^Zi and xFsz, in lowest order. In 0{as) also the longitudinal 
structure functions 



Sl{x, Q') = S2{x, Q') - 2xSi{x, Q'), 



(3) 
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with S = F,G, H,W^ , and Fz, contribute to the scattering cross sections. 

In the case of polarized lepton-polarized hadron scattering a similar amount of 
structure functions occurs ^. In the kinematical range of the present experiments 
the scattering cross section is determined by the I7P term, however, and in this 
approximation only the structure functions gi[x,Q'^) and g2{x,Q^) contribute. 

Not all of the structure functions mentioned above can be determined at a suf- 
ficient accuracy combining different cross section measurements ^. The structure 
functions F\'^''^'^\x,Q^), and the combinations W^^{x,Q^) = | (^"2^ + ^"2") , and 

xW^^^x^Q"^) = I (a^VFs'' + a^VFg"^ can be measured precisely in two variables [x^Q"^) 
in a wide kinematical range. For them QCD analyses may be performed to determine 
the QCD parameter A and constraints on the gluon density. 

The measurement of the longitudinal structure functions Fj^ ^^'^\x,Q'^) is more 
difficult and requires a variation of the CMS energy keeping x and fixed. A 
precise determination of this structure function is of special importance since it is 
directly related to the gluon density, see sect. 4.2. For other structure functions as 
xG^2 p(<i)^ i^2Z5 only the a;-shape may be determined since they are measured from 
cross section differences, or, the reconstruction of the kinematical variables turns out 
to be difficult. 

The determination of the QCD and QED corrections to the different deep inelastic 
scattering processes is of great importance for the quantitative understanding of the 
nucleon structure. In the present paper we give a survey on the status of perturba- 
tive QED and QCD radiative corrections to deep inelastic scattering processes. We 
also include a discussion of the status of the determination of a^, and of the QCD 
corrections to some exclusive processes, as the heavy flavour structure functions and 
J/^ production in IN scattering, which play an important role for the determination 
of the gluon density of the proton. 

2. QED Radiative Corrections to DIS 

The QED radiative corrections to the deep inelastic scattering cross sections may 
become rather larger in some kinematical ranges. They have to be precisely known 
to unfold the neutral and charged current structure functions from the scattering 
cross sections. The first dedicated calculation of the radiative corrections to deep 
inelastic eN scattering was performed by Mo and Tsai ^ and used in the analysis 
of the SLAC experiments. Later calculations were performed in refs. for l^N 
scattering. The detailed knowledge of QED and electroweak radiative corrections was 
of special importance also for the measurements of the electroweak parameters in deep 
inelastic v{y)N scattering With the advent of HERA the radiative corrections 

were partly recalculated and dedicated calculations for deep inelastic neutral and 
charged current e^p scattering were carried out by different groups using different 
techniques These approaches include both semi- analytical calculations ^^"^t 
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and calculations based on Monte Carlo techniques ^8-30 

Dominant contributions to the QED radiative corrections may be obtained using 
leading log (LLA) techniques i^~24_ This approach, which is based on the factorization 
of (coUinear) fermion mass singularities, allows to determine the terms oc a\n[Q'^ / m"^^) 
in a straightforward way for different settings of the measured kinematical variables. 
Also higher order terms were calculated within this approach 21,24 fpj^g LLA QED 
radiative corrections may be described by 

dxdy dxdy dxdy dxdy dxdy dxdy dxdy 

Here, I labels the type of bremsstarhlung, which is in lowest order either initial or 
final state radiation. In higher orders also products of initial and final state radiation 
terms contribute. In many situations final state radiation does not occur, or can be 
delt with in a cumulative way due to the calorimetric measurement of the scattered 
electron. Then all terms refer to initial state radiation only. The label C denotes 
the Compton-contribution" i.e. the coUinear term formed by low radiation of the 
virtual photon from the initial state hadron or quark lines. The soft contributions, 
> 2, soft, can be exponentiated. Starting with (9(a^) also terms due to e~ — )• e"*" 
conversion are present in the leading logarithmic order. 
The first order terms are described by: 



^ie/ dzP^ine{z-z,)J{x,y,Q') 



dxdy 27r y ee dxdy 

where 



dxdy j ' 
(5) 



P^l' = ^ (6) 
i — z 

denotes the non-singlet QED splitting function of a massless fermion into a fermion. 
The scale of the correction is set by the logarithm 

Le = In^ - 1. (7) 

This notion reproduces the soft photon terms of complete calculations in leptonic 

variables (cf. e.g. ^^). The shifted variables and the Jacobian J7 depend on the 

choice of the outer kinematical variables (see for a summary of these terms). 
The Compton term is given by 



d 



a' 



dxidyi xiS 



(8) 



"This term was already found in ref. 
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for leptonic variables in LLA. A more refined expression was derived in ref. LLA 
second order corrections are easily obtained by convoluting with the leading order NS- 
splitting function. Although the Compton-type contribution counts to the radiative 
corrections to deep inelastic scattering in an inclusive description, its experimen- 
tal signature is rather different compared to typical deep inelastic events, showing 
a photon-electron pair which is nearly balanced in p± and little hadronic activity 
only Such a signature can be easily tagged. Due to this one may even use these 
events to measure nucleon structure functions both at small x and small Q^. 
The second order corrections (9((aie)^) are: 



J2 ^l,2loop 

dxdy 



27r 



jyzP2''\z)\e{z-Zo)J{x,y,z) 



d^a° 



dxdy 



dxdy 



',S=S 



+ (f dz {lIp(^'%) + ie E In ^ pS\z)] 
\ZTry Jzo f^i^g rrif J 

Here the different second order splitting functions are given by 







dxdy 



hs=s 



(9) 



2 

l + z2 r 



21n(l - z) - Inz + 



Pt'\z) = ^ [i^,^) ® P«] (.) 



+ -(l + z)lnz-(l-z). 



{l + z)\nz + ^{l-z) + jkl-z% 



pS\z) = N.if)Q}lpil\z)e{l-z-'-^) 



(10) 

(11) 
(12) 



denoting double-photon radiation, scattering of a fermion into a fermion by a coUinear 
photon, and coUinear fermion pair production. (/) denotes the Mellin convolution 



A{x) (x) B{x) = I dxi I dx28{x — xiX2)A{xi)B{x2) , 
Jo Jo 



(13) 



Qf IS the fermion charge, and Nc{f) = 3 for quarks, Nc{f) = 1 for leptons, respec- 
tively. 

The soft-photon exponentiation is performed solving the non-singlet evolution 
equation in the range z — )• 1 analytically (cf. e.g. ^^). Since the terms up to (9(a^) 
were taken into account in eq. (9) already the corresponding contributions have to be 
subtracted. One obtains 



£^{>2,soft) 

dxdy 



dzP^>'\z, Q') I e{z - Zo)J{x,y, z) 



dxdy 



,s=s 



dxdy 



(14) 



5 



with 



Ztt l — z k Ztt 
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+ 21n(l - z) 



L>Ns{z,Q ) - C(l-^j _^ 5 

C = -31n[l-(a/37r)ie]. 
Finally, the fermion conversion term in 0{q^L\) reads 

^2^(2,e-^e+) i ^2^(0) 



/" dzP{z,Q'^;e e'^)J{x,y, z) 

J Zn 



dxdy 



dxdy J zo 

with the conversion rate given by 

P{z,Q^;e-^e^)=[^yLlp(^'%) 



(15) 

(16) 
(17) 

, (18) 



,s=s 



(19) 



Lepton Measurement 
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Figure 2: a) Oipt) initial state QED corrections to neutral current deep inelastic scattering for 
leptonic variables in LLA. b) initial state QED corrections using the double angle method. A 
z-cut o\ E = 35 GeV (cf. ^^) is applied. Full lines: LLA corrections up to second order + soft 
exponentiation; dashed lines: 1st order LLA corrections. 

In figure 2 we illustrate the size of QED radiative corrections for two types of the 
measurement of the kinematical variables in the case of neutral current deep inelastic 
scattering. In the case of the (classical) leptonic variables x and are determined 
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from the measurement of the scattered lepton only. In the double angle method the 
scattering angles of the outgoing lepton and hadronic jet are used on the other hand^ 
Whereas the radiative corrections for leptonic variables become very large at small 
X and high y, they behave flat in y in the case of the double angle method and are 
very small in the range x < 0.01, where the structure functions become rather large. 
Thus the latter set of variables behaves ideal in this respect. 

The different methods to calculate the QED radiative corrections have been com- 
pared in 0{a) for a variety of kinematical measurements and are well understood. In 
2nd order so far only LLA results are available ^^'^^ for the full set of outer kinematical 
variables studied by the HERA experiments. 



3. The running coupling constant 



The strong coupling as{fi'^) is a central parameter in QCD. It is not an observable 
itself but the various hard scattering processes are often compared with respect to 
this quantity. Since it is scheme-dependent these comparisons have to be performed 
in a single renormalization scheme, as e.g. the MS-scheme The running of as{fi'^) 
is determined by the renormalization group equation 



/3o 2 



(20) 



ln/i2 47r ' (47r)2 ' [ATrf 

So far the contributions to the /3-function have been calculated up to 3-loop order in 
the MS-scheme, where the LO NLO and NNLO terms are given by 

2 



/3o 

^2 



11 - 








102 - 


- Y^/, 


2857 


5033 


2 


18 



325 , 

Nf + 



54 " ^' 

Nf denotes the number of active flavours. The solution of (20) reads: 

^r^)^^Am)^ A) - .^")(«.(g^); A 



(21) 
(22) 
(23) 

) (24) 



The superscript n denotes the term at which the expansion of the /3-function in (20) 
was truncated. In NNLO one obtains 



87r/3o 



In 



Wtt^x 



2^2 



167r2/3o + 4/3i7ra; + (32X^ 



PI - 2^2 



arctan 



27r/3i + 



87r/3o74/32/3o - P! VS^VwT^ 



(25) 



For other types of measurements see and references therein. 
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Note that 



Nf<5: 4/3o/32 - Z?' > 
Nf = 6: 4/3o/32 - f^l < 0. 



(26) 



Table 1: 

Summary of recent measurements of (from ^^) and results from lattice calculations (cf. ^^). 





(Q) 




a 


s{Q) 




Process 


[GeV] 


«.(M^) 


exp. 


theor. 


Theory 


GLS (CCFR) 

(CLEO) 
R^ (ALEPH) 
Rr (OPAL) 
R-r (P. Raczka) 
7?e 77 (CLEO) 

T(15) (CLEO) 
ep jets (HI) 
PP Wjets (DO) 


1.73 

1.78 
1.78 
1.78 
1.78 
2.98 
9.46 
5-60 
80.6 


+ 0.007 
- 0.009 
0.116 ± 0.003 
0.122 ± 0.003 
0.123 ± 0.003 
0.120 ± 0.003 
0.101 ± 0.010 
0.111 ± 0.006 
0.123 ± 0.018 
0.121 ± 0.014 


+ 0.006 
- 0.007 

0.002 

0.002 

0.002 

0.002 

0.008 

0.001 

0.014 

0.012 


+ 0.004 
- 0.006 

0.002 

0.002 

0.002 

0.002 

0.006 

0.006 

0.010 

0.005 


NNLO 

NNLO 
NNLO 
NNLO 
NNLO 
NLO 
NLO 
NLO 
NLO 


DIS: 

pF2,xFs 

flF2 


5 

7.1 


0.111 ± 0.006 
0.113 ± 0.005 


0.004 
0.003 


0.004 
0.004 


NLO 
NLO 


e+e- ^ jets (CLEO) 


10.53 


0.113 ± 0.006 


0.002 


0.006 


NLO 


e+e- ^ 

scaL vioL (ALEPH) 
ev. shapes (SLD) 


91.2 
91.2 


0.127 ± 0.011 
0.120 ± 0.008 


0.003 


0.008 


NLO 
resum. 


7(Z° ^ had.) (LEP) 


91.2 


0.127 ± 0.006 


0.005 


+0.003 
-0.004 


NNLO 


LOT 




0.111 to 0.115 




0.005 to 0.007 


43 



One may solve eq. (20) using as{Qo) = Ois{Mz) as input. Due to the fact that Nf = 5 
for Q < 10 GeV the NNLO correction diminishes the NLO solution slightly at low 
values of cf. 

In table 1 recent measurements of (cf. ^^) from different Q ranges are compared, 
which were evaluated at the Z° scale. We also added the range of recent results from 
lattice calculations (LGT) (cf. ^^). The most precise measurements stem from the 
various high statistics deep inelastic scattering experiments at the one side, and high 
precision measurements of different observables in e"'"e~ annihilation, on the other 
side. The average values for these measurements are: 

DIS : a;{Mz) = 0.112 ± 0.004 (27) 
e+e" : a:{Mz) = 0.121 ± 0.004, (28) 
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showing a 2<t difference at present. The results from lattice calculations yield val- 
ues of as{Mz) = 0.111. ..0.115. Here, the systematical error is estimated to be of 
(9(0.005. ..0.007) still Whereas the pure gauge-field contributions are rather well 
understood the quark terms deserve more detailed investigations in the future to 
obtain decisive results. 



4. The Evolution Equation 



The structure functions, Fj[x,Q'^), describing the deep inelastic scattering cross sec- 
tions may be expressed as a convolution of the bare parton densities, and hard 
scattering cross sections, <t^-. 



Fj{x,Q^) = fi{x)(g)ai{a,,Qyn\x,e). 



(29) 



The functions a'j contain initial state mass- and ultraviolet singularities. If the calcula- 
tion is performed in 4 — e dimensions they emerge as poles in e. The mass singularities 
can be factorized and absorbed into the bare parton densities The ultraviolet sin- 
gularities are removed by the renormalization of the bare coupling constant, a^- In 
this way two scales, the mass factorization scale, /ii, and the renormalization scale, 
/i2, are introduced. One further separates ctj into its pole-, F^, and non-pole parts, 
Cjk- This separation is obviously arbitrary and introduces a scheme-dependence. One 
obtains: 

Fj{x,Q') = f\x)<g)T^{a,{^il),^il/^i\^il/^il,e)<g)CJkM^^^^^ (30) 

The first two factors on the rhs of (30) define the renormalized parton densities, 
which are scheme-dependent. One may identify the scales fii = fi2 = M and obtains 



Fj{x, Q') = fix, a,{M'), M'/fi') Cjk{as{M'), Q'/M\x). 
Let us transform the above equations to moment space by 



r dxx^-^A{x) = A{N). 
Jo 



(31) 



(32) 



The invariance of the structure functions Fi against the choice of M may be expressed 

by 

^5^+^^^^^"^^'^^^^ ^.W = (33) 



leading to 



M^+f3(g)^--y^(g) 



fk{N) 
CAN) 



0, 
0, 



(34) 
(35) 
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where we use a generic notation for the anomalous dimensions, 7^. Here, eq. (34) 
describes the evolution of the parton densities, which results directly from the renor- 
malization group equation. The anomalous dimensions 7^ are the Mellin transforms 
of the splitting functions, which will be defined below. 

Let us introduce some combinations of parton densities: 



q^{x,Q') = Y:qt{x,Q') 



(36) 
(37) 

(38) 
(39) 



ql and q^ denote flavour non-singlet combinations, and q'^ is the singlet density. 
The Mellin transform of (34) yields the non-singlet and singlet evolution equations in 
X space: 



<1^{^,Q') = ^4^P-{x,a,)^q-{x,Q') 



d 



dlnQ' 



d 



dlnQ^ 



q^x,Q') 
G{x,Q') 



2tv 
27r 
27r 



P+{x,as) (g) qi {x,Q^) 



P{x,as) ® 



q^x.Q') 
G{x,Q') 



(40) 
(41) 
(42) 



The splitting functions, and P can be expressed by the following perturbative 
series: 



PHx,a.) = Pffi + gP"'«(-)+(g)'p"'(^)(-) + 
P(^,a.) = p(°) + ^pW(^)+f^)'p(^)(^) + ... 
Let us change the evolution scale Q = M introducing 

t 



-dlnQ' 



+ 



dt. 



27r " ~" \ 2/3o 27r 

Then, the evolution equations may be rewritten as: 

q^{x,t) := E-{x,t)(g)qr{x) 

qt{x,t) := E+{x,t)(g>qt{x,t) + ^[Err{x,t)-E+{x,t)](g>q+{x) 



(43) 
(44) 

(45) 
(46) 

(47) 
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E{x,t) (g) 



(48) 
(49) 



where we introduced the evolution operators and E obeying the initial conditions 



limE^(a;,^) = ^(1 - a;) 



YimE^{x,t) = 1^(1 -a;). 



(50) 
(51) 



They allow to separate the non-perturbative input densities qt{x),q'^{x) and G{x) 
from those terms which can be calculated perturbatively. The evolution operators 
itself obey the evolution equations: 



j^E^M = Ipns{x) + ^R^{^) + •••} ® 



d 
di 



E{x,t) = !^P^°\x) + ^R{x) + ...| ® E{x,t), 



where 



R^{x) 



R(x) = pi^)(x) - ^P('^ 



(52) 
(53) 

(54) 
(55) 



^.1. Splitting Functions 

The well-known leading order singlet splitting functions are given by^^~^° 



— + -8{l-z) 

(l-z)+ 2 



i + {i-zy 



gg \ ! 



1C, 



z 

1 - z 



+ 



(1--)h 



+ -/3o^(l - z\ 



(56) 
(57) 
(58) 

(59) 



and the non-singlet splitting funtion obeys P^\{^z^ = P^^\z). Here, Ca = Nc = 3, 
Cp = {N^ — l)/iVc, and Tf = 1/2. In leading order the splitting functions for space 
and timelike virtualities are the same This relation is violated in higher orders. 
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The non-singlet splitting functions in NLO for spacelike virtualities are given 

by: 

P±(z,a,) = P^{z,as)-S{l- z) f' dzp-{z,as), (60) 

P'^{z,as) = Pqq{z,as) ± Pqq{z,as), (61) 

where the superscript it labels the type of non-singlet evolution, see eqs. (40, 41). 



+ [cIPf{z) + ^CfCaPg{z) + GFNfTfP^^iz) 

Pa{z) 



(62) 
(63) 



1 -l- z'^ / 3 \ 1 

Pf{z) = -2— — lnzln(l-z)-l + 2z I In z - -(1 + z) W 



Pg{z) 
Pn,{z) 

Pa{z) 



-5(1 -z) 
l + z2 



(64) 
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In z H In z + 



67 



-TT 



40 



+ 2(l + z)lnz + — (1-z) (65) 



i±^(lnz + ^)+2(l-z)j (66) 
2-^/ —In +2 1 + z lnz + 4 1 -z 67 

1-z Jzl{l+z) U \ U / 



The NLO singlet splitting functions were derived in refs. for the unpolarized 

case both for space- and timelike virtualities. 

^.2. Coefficient Functions 

In the calculation of the higher order corrections to the structure functions the coef- 
ficient functions are required (cf. eq. (31)). These are scheme-dependent quantities. 
In the MS-scheme they read in O^a^) 



c?Az) 

G?Az) 
C«(z) 

G^gKz) 

G^gKz) 



1 + z^ / 1 - -2 3\ 1, 



c?Az) - Cf{i + z) 

2NfTf I [z' + (1 - zf] In - 1 + 8z(l - z)| 

cSkz) - 8NfTfz{l - z). 



(68) 

(69) 
(70) 

(71) 
(72) 
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The 0{al) contributions to the coefficient functions were calculated in refs. 
for the structure functions F2, Fl, and XF3. Whereas in refs. the coefficient 

functions were derived in z-space, the moments M"|^_2 ^^re calculated for Fl 
and F2 in ref. The results of both calculations do fully agree. 

To illustrate the numerical importance of the 0{al) calculation we compare in 
figure 3 the 0{as) and the 0{al) result for the structure function Fl[x,Q'^) in the 
MS-scheme using the parametrization for the parton densities. 



^ 0.8 

!-f 0.7 





Q' = 20 GeV^ 






van Neerven, MSbar 






moss factorlzotion, /i^ = (f 






LO, CTEQLO 






NTLO, CTEQ2M 
1 1 j~r-r 






10 10 10 ' 


X 



Figure 3: Comparison of QCD corrections to the structure function Fii^x^Q^). 0{as) (dashed 
line); 0{al) (full line); from 



The NLO correction leads to a depletion relatively to the LO result by ~ 15% at 
X ~ 10~^. In the small x range Fl[x,Q'^) is widely determined by the gluon density. 
The 0{a.l) corrections are not negligible and have to be taken into account for the 
unfolding of the gluon density from a measurement of Fl[x, Q^). 



4.3. 0{al) corrections 

For a series of observables related to deep inelastic scattering QCD corrections up to 
3-loop order were calculated. The quantities which were studied so far are moments 
of structure functions and specific combinations which are related to sum-rules 
of structure functions. 

The 0{a.^) corrections to the Bjorken-^^, the Gross-Llewellyn Smith-™, and the 
polarized Bjorken^^ sum rules are: 

l^x [Fr{x,Q')-Fr{x,Q')] = 1-^^-2.3519 (^)'-8.4852 (^y + ... (73) 
l'dx[FrM^)-FrM^)] = 6{l-^+(^)'(-g + ^iV,) 
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+ 
+ 

+ 

Note the small difference in the correction factors between the Gross-Llewellyn-Smith 
and the polarized Bjorken sum rules in the Nf term in 0{a.^). 

These calculations are performed by fast formula manipulation programs as FORM 
Due to the complexity of the problem they request hundred(s) of CPU hours on 
present day computers. The non-singlet moments M"|„=2,4,6,8 of the structure func- 
tions F2{x,Q'^) and Fl{x,Q'^) were also calculated The calculation of still higher 
moments for the non-singlet case and the first singlet moments of these structure 
functions is being performed currently 

Precise experimental data for the sum-rules and moments quoted will allow very 
concise tests of QCD in its perturbative range. 

5. Resummation of small x contributions 

At small values of x contributions to the splitting functions show a singular behaviour. 
This is both the case for the singlet and non-singlet splitting functions. Resummations 
of large terms are also required at large x Here we will deal with the small x 
behaviour only. 

5.1. Singlet terms 

In the singlet case the singularity in the splitting function is of the type P oc 
[l/x)a[ In ~^ X. In leading order these contributions are resummed by the BFKL equa- 
tion The leading singular terms in the gluon anomalous dimension, 7^^, to all 
orders in are given by the solution of 

N-1 = a,x[7L(^,«.)], (76) 
X{z) = 2^(l)-^(z)-^(l-z), (77) 

with CKs = CaO-sIt^- The solution of (76) is multivalued and one has to select a single 
Riemann-sheet. This is done imposing the condition 

lim 7i(iV,a7)oc-^, NeC. (78) 

|iV|^oo iV — 1 

"^For a recent review see ref. "^"^ . 
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216 

, 10009 91 , 

Nf h — Cs 

^ ^ 1296 54^ 



44 55 
y Cs + yCs 



9a 



9v 



1--. 
L TT 



■Nf 



10339 61 , 

\ Cs 

1296 54^ 



(75) 
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In figure 4 real and imaginary part of this solution are shown. Asymptotically, i.e. 
for A = as/{N — 1) ^ 1, one obtains: 



7L(iV, as) 



a. 



N 



A + 2CsA^ + 2CsA^ + 12CIA' + 



(79) 



The solution of (76) using (78) contains three branch points. They are solutions of 
the equation (z = 7^): 

Z Sm TTZ 



and are ^^'^^ given by: 



ill = 1/2 
7l2 3 = -0.4252 ± 0.4739i 



(81) 
(82) 



The former value is seen in figure 4 as the 'roof in Re7, the latter ones are the edges 
at p = -1.4105 ± 1.9721i. 




S ' 

0.75 
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0.25 
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-0.5 
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Figure 4: The solution of eq. (76) for complex values of /9 = (iV — l)/as. a) Re(7i); b) lm(7i), cf. 

The resummation for the most singular terms in the singlet anomalous dimension 
takes the form *° : 

- £ (^)'4? + t^. + O U {^y) . (83) 



79 



The LO and NLO terms to this matrix equation are given by 

/ \ 

r^(^) = [ §j7L(iV) MN) ) 



(84) 
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(85) 



So far only the quark contributions in the NLO terms are calculated. These, and the 
LO terms, are functions of 7^. An approximate solution for 7jvl reads ^""^ 



2a. „ [ ocs I OL 

jnl{N) ~ -^Tf \ 1 + 2.17^^ + 2.30 



37r 



N - 1 



N - 1 



+ 8.27 



a. 



N -1 



+ 



(86) 



One may use the expression (83), supplemented by the terms contributing to {"fab) 
up to NLO (cf. section 4), and solve the singlet evolution equation. This was done 
in ^^'"^^ assuming a flat input at Ql = 4GeV^. 





HI 1993 data 
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Figure 5: Resummed predictions for the structure function F2{x,Q'^) see text. 

In figure 5 the behaviour of F2{x,Q^) due to different contributions to the evolution 
is illustrated. The evolution of the flat input using the pure LO and NLO terms leads 
to a small rise at low x. The resummed LO terms affect F2 through the gluon density 
only and lead to a correction which is rather small if compared to the correction due to 
the resummed NLO terms, which affect the quark densities directly. Here, for the yet 
unknown entries 7^^^ in eq. (85) the NLO terms were used. Since the resummation (83) 
is performed for the most singular terms at small x only, momentum conservation has 
to be restored explicitely. This is possible in different ways. One may either introduce 
an appropriate term oc ^(1 — z) in the splitting functions or multiply the anomalous 
dimensions by (1 — N). The latter solution is illustrated by the full line in figure 5, 
while the former is shown as the dotted line. The large difference between the two 
solutions shows that it is likely that sub-dominant terms may be as important for the 
scaling violations of F2 as the singular terms at small x. 

"^Related numerical studies were performed in ref. ^^'^^ ^Iso. 
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A unified form for a leading order evolution eqution accounting both for the small x 
terms due to the BFKL equation and the evolution kernels was found in 

It is based on the angular ordering of the gluon cascade. In the limit of small x the 
BFKL equation is obtained, while for medium and large values of x this equation 
turns into the LO evolution equation of pure gluodynamics, i.e. the angular ordering 
turns into strong ordering of ki_. This concept was worked out in leading order so 
far. A numerical illustration of this resummation was given in 



5.2. Non-singlet terms 



As shown in sect. 4.1 the most singular terms in the non-singlet splitting functions 
in 0{as) and 0{a1) behave as as{as\rL' xY . In ref. a resummation of these con- 
tributions to the lb combinations (cf. eq. (60)) of non-singlet structure functions was 
derived. Recently very sizeable corrections due to this resummation have been 
claimed both for unpolarized and polarized structure functions. 

Similar to the considerations in section 5.1 the resummation can be studied in 
the context of the renormalization group equation, see ref. Unlike the case of the 
BFKL equation the present resummation deals not with the anomalous dimension 
but with the structure function itself. Therefore one has to consider the evolution 
equation for the non-singlet structure functions here: 



where a. 



das poai 
as(Q^)/(47r). In NLO the evolution kernels K^g-^ are 



Po 



(87) 



(88) 



771 en 

^2 5 



The labels it denote the type of the non-singlet evolution. The combination 
e.g., belongs to the '+' type, and xF^^ -\-xF^^ and g^i —g^ are '— ' type combinations. 
In the latter case the splitting functions obey Jq dxPfi^x) = in each order in a, due 
to fermion number conservation. 

The resummation of the most singular parts in the kernels Kf{^x.,aa) read in 
iV-space 



-2N I 1 



-2N I 1 



2a,Cf 



^ ttN^ j 



(89) 



\ 



ttN^ 



2Ncas d 
ttN dN 



In (e^'/41)-i/2iV|(^)) 



(90) 



where N = z<J'a^/2, and = Ncas/ir. Dp[x) denotes the function of the parabolic 
cylinder. Up to 0{al) the coefficient functions in the MS-scheme behave at most 
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oc In^ ~^ X. Therefore the contributions to Kj^g ^^^i^x^as) can be directly compared 
with the results from fixed order perturbation theory at least up to NNLO in the 
MS-scheme. These are known up to NLO and are found to agree 




Figure 6: The small-a; evolution of non-singlet structure functions in NLO and the correction 
due to the resummed kernels eqs. (89,90)^°: a) F^^ - b) gl^ - gl'^. The labels A 

and B refer to the way in which fermion number conservation is restored. A: r~(iV, a^) 
r-(iV,a,)-r-(l,a,), B: T- {N, a^) ^ T- {N, a^) ■ {I - N). 

In figure 6 we compare the evolution for + and — non-singlet combinations of structure 
functions in NLO with the resummed contributions beyond NLO. The latter terms 
yield corrections at the level of (9(1%) in the accessible kinematical ranges (cf. ^^). 
Moreover, in the case of — type combinations different possible ways to restore fermion 
number conservation lead to a variation of the terms beyond NLO by a factor of three. 
This signals that yet unknown medium- a; contributions to the evolution kernels may 
be as important as the small x terms. A similar behaviour was observed in the case 
of the singlet terms in sect. 5.1. 



6. Heavy flavour contributions to structure functions 



In lowest order the contributions to the heavy flavour structure functions F2,l{x, Q^) 
for the term in neutral current deep inelastic scattering are described by the 
diagrams due to photon-gluon fusion. They are given by 



F?2{x,Q\ml) = 2elx^^ f 

I'K Ja 



ax y 



where 



9,2 



Ml 



1 

2 

+ /3 



+ {1- zf + z{l-3z) 



X mg 



(91) 



-1 + 8z(l -z)- z{l - z)- 



4m2 



In 



1+/3 
1-/3 



(92) 
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Ami 1+B 



/3 



(93) 



with a = 1 + AM^/Q"^ and /^^ = 1 _ (4771^(5^)2(1 - z)-^ These contributions and 
the other LO terms were derived in refs. ^^-99 

The NLO contributions were calculated in ref. loo-io^ pjienomenological studies 
can be found in i°4,i05 NLO corrections stabilize the numerical values of the 

scattering cross section with respect to the choice of the factorization scale. 

In figure 7 numerical results are shown for the LO and NLO contributions to 
F^'^i^x^Q^) illustrating the size of the NLO corrections vs the LO term. 





Figure 7: The a;-dependence of i^^'^''^'^(a;, Q^, ttz^) (lower pair) and F^^^''^'^(x,Q'^,ml) (upper 
pair) at fixed Q^. The solid lines are for = 100 GeV^ and the dashed lines are for = 
lOGeV^; (from i°i). 

7. J/ip production 

As in the case of heavy flavour production from the measurement of the deep inelastic 
production cross section of J/^ particles the gluon density can be determined. In 
lowest order the cross section for the photoproduction case is given by 

da^^ 128^0^ \m\'s'sl + tHl + u'ul 

in the colour- singlet model. Here we used the abbrevation Vi = r — Mj^^^r = s,t,u. 
The scattering cross section for finite photon virtualities [Q"^ > 0) was derived 
in At a photon beam energy of Ej = 150 GeV, e.g., the production cross 

section is lowered by a factor of ~ 7 for = 20 GeV^ compared to = 0. 

Recently the NLO corrections to the photoproduction cross section have been 
calculated i°9,iio_ stabilization of the scale behaviour in the range Q'^/ml > 1.5 
was obtained in comparison with the leading order result. 
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Figure 8: The total cross section -\-p J/^) + X) as a function of the photon-proton CMS 
energy for different parametrizations of the parton densities (from ^^°). 



In figure 8 the total production cross section (NLO) is shown and compared with 
recent measurements at HERA. 



8. QCD corrections to polarized structure functions 



mg 



The leading order singlet splitting functions for polarized deep inelastic scatter- 

111,112,48 



are 





— Cp - 








— Cp 


p^'M 


= - 



1 + 3 
. , +-8il-z) 



(1 - -f] 



z 
2 



4z + 2 



(95) 
(96) 
(97) 
(98) 



Again the non singlet splitting function obeys Pns{x) = Pqq{x). The non-singlet 
splitting functions in NLO are known from the unpolarized case (see sect. 4.1) al- 
ready. Recently also the singlet splitting functions in NLO were calculated in the 
MS-scheme. The result of this calculation has been confirmed in recently. It is 
interesting to note that eqs. (96, 97) differ from eqs. (57, 58) by relative signs in some 
of the terms. A common characteristics of these quantities is their leading singularity 
behaviour at small x. Unlike the unpolarized case the Mellin transform of the leading 
terms behave ~ 

The LO coefficient functions in the MS-scheme are^^^~^^^ 

(ln{l-z)\ 



a. 



8[l -z) + :^Cp 



1 + 

2^- ln(z) + 4 + 2z - ^(1 - z)(4C(2) + 9) 



2(1 + z)ln(l - z) 



(99) 
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C,{z) = -^NfTf {4(2z - 1) [ln(l - z) - ln(z)] + 4(3 - 4z)} . (100) 

The NLO coefficient functions were calculated in ref. 

With these quantities at hand the scaling violations of the structure function 
g\{x,Q^) can be studied up to NLO. This has been done recently^^^'^^° extending 
earlier studies in leading order. Still the parametrizations of the polarized parton 
densities do widely vary. Particularly this holds for extrapolations to the small x and 
large ranges. 




Figure 9: LO Parametrizations of the structure function gi{x,Q^) in the range a; > 10 ^. Full 
line: = IQ GeV^ dashed line: = GeV^ dotted line: = lo^ GeV^ dash-dotted 
line: Q"^ = 10^ GeV^ (cf. ^^); a) parametrization ref. b) parametrization ref. 

A systematic collection of the parametrizations for polarized parton densities and 
their evolution in LO (later than 1989) can be found in 

9. Open Problems 

At the end of this survey I would like to list a series of open problems, the solution 
of which is of importance for the forthcoming development of the understanding of 
deep inelastic scattering in the perturbative range. 

Concerning the QED corrections Oi^a^Le) terms should be calculated for sets of 
variables in which the 2nd order LLA terms are still large. Furthermore, complete 
(9(a) calculations should be performed for some sets of kinematical variables for which 
only LLA results exist to further improve the numerical accuracy. 

The knowledge of the 3-loop non-singlet and singlet splitting functions for un- 
polarized deep inelastic scattering would be important to perform QCD tests at the 
level of NNLO corrections using measurements of the structure functions F2{x,Q^), 
^^2(3^5^^)5 and xW:i{x,Q^). As a by-product of these calculations also a partial test 
of the validity of the resummations discussed in section 5 would be obtained. 

The calculation of the NLO corrections to the BFKL equation are needed for 
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futher studies of the behaviour of the anomalous dimension at small x. 

NLO corrections should be calculated for different processes in polarized lepton- 
polarized nucleon scattering. They are also not yet performed for deep inelastic 
leptoproduction of J/^ particles at > 0. 

Besides of the twist-2 contributions to structure functions being discussed in the 
present paper the understanding of higher twist terms is important. Particularly 
gluonic twist-4 contributions to F2{x,Q^) would be interesting to be derived in a 
complete calculation to compare with results obtained in Regge approaches for the 
limit of small x. These terms may play an important role in understanding the 
screening of the gluon density. 
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